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Abstract 

Tunnel splitting oscillations in magnetic molecules are reconsidered within the simplest model 
for the problem, which does not contain fourth order anisotropy. It is shown that at large magnetic 
field, there is only one instanton, and it is continuous. This is in contrast to the discontinuous 
instantons that are induced by the fourth order term [Ersin Kegecioglu and A. Garg, Phys. Rev. 
Lett. 88, 237205 (2002)]. 
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I. INTRODUCTION 



The purpose of this short note is to clarify the differences between certain types of instan- 
tons that arise in the study of spin orientation tunneling in molecular magnets. The problem 
derives from the specific phenomenon of tunnel splitting or gap oscillations in the molecu- 
lar ion [(tacn) 6 Fe802(OH)i 2 ] 8+ (abbreviated to Feg henceforth) [[[[], but our considerations 
should apply to other molecular magnetic system as well. 

The simplest model Hamiltonian that leads to gap oscillations is 

n = k,j 2 z + k 2 .P y - gf i B j ■ h, (i.i) 

where J is a dimensionless spin operator, k\ > k 2 > 0, and H is an external magnetic field. 
It is known that for H||z, the tunnel splitting, or gap, between the two lowest eigen states 
of 7io oscillates as a function of H z , vanishing whenever || 

H z = j(j-n-f)H c . (1.2) 

Here, 

A = k 2 /ki, and H c = 2kiJ/g^B- (1-3) 

For Fes, for which J = 10, k\ = 0.338 K, and k 2 = 0.246 K, this would imply that the 
splitting is quenched at 10 values of H z > 0, equally spaced by 0.263 T. This quenching 
phenomenon is nicely interpreted in terms of interfering semiclassical paths, or instantons 
I- 

In reality, only four quenches are seen in Feg, and the spacing between the quenches is 
~ 0.41 T, 50% larger than predicted by Eq. ( |1.2|) . Both these facts are accounted for by a 
slightly modified Hamiltonian 

H^Ho- C[(J Z + iJ y ) 4 + h.c], (1.4) 

with C = 0.29 /iK, as is easily verified by explicit numerical diagonalization. The fact that 
only 4 quenches are seen is paradoxical, however. The relative phase between the interfering 
instantons is a topological quantity, so that it must equal J-zr at H = 0, and tend to as 
H — > H c . There is no apparent way to get other than 10 quenches. The resolution is that 
the C term, though small, has a dramatic effect: any non zero C leads to the existence 
of instantons that are discontinuous at the end points [§J. One of these instantons ends 

2 



up being the dominant one by virtue of having the lowest action beyond a certain field 
(~ 0.25if c for Fes), and since it has no interfering partner, there are no more quenches of 
the splitting beyond this field. In this way, the discontinuous or jump instantons may be 
said to shunt the interference effect. The jump instantons have discontinuities, and thus do 
not have to respect the topology. 

In private correspondence to us, E. Chudnovsky has alleged that jump instantons were 
first introduced by him and Hidalgo in an earlier study of the simpler Hamiltonian ( |1 . 1|) 
at high fields ||. We disagree. The instanton studied by Chudnovsky and Hidalgo (CH) 
is of exactly the standard continuous type arising in the study of a particle tunneling in a 
symmetric one- dimensional double well, but its nature is obscured by the subtleties of spin 
trajectories, and the singularities of the spherical polar coordinate system on the sphere. 
These subtleties are vexing and the source of innumerable pitfalls, so it seems useful to us 
to elaborate on the instantons in question, and show explicitly that they are continuous. 
We also use the opportunity to correct a misconception in a still earlier study of the same 



problem, i.e., Eq. ( |1.1| ) at high fields, by one of us 



II. DETAILS 

We consider only fields parallel to z, and define the reduced field 

h = H/H c , (2.1) 

and the special value 

h* = (1 - A) 1/2 . (2.2) 

CH study the problem for h > h*, and Chudnovsky is asserting that the instantons in that 
study have jumps. We give three arguments to show that this is not so. 

The first and mathematically simplest way is to consider fields just less than H c . We 
write 

h = (l-v) (2-3) 

where < r\ <C 1. The expectation value of J z in the low energy states is then extremely 
close to J, so we can approximate 

J x ~ Jq, J y ~ P, (2.4) 
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where p and q are momentum and position operators with the commutator [q,p] = i. It is 
self-consistently verifiable that J y ~ rj, J x ~ r/ 1//2 , so that we can write J z ~ (J 2 — J 2 — J 2 ) 1 / 2 . 
Expanding the term —gfj, b HJ z to order J 2 and J*, rewriting in terms of q and p, and omitting 
an additive constant, we arrive at the equivalent Hamiltonian 

H equiv = ^p 2 + h.hJY - rjhjy. (2.5) 

Here, 

m= [2(k 2 -ijk 1 )] - 1 « (2A&)" 1 . (2.6) 

Equation ( |2.5| ) is the Hamiltonian for the particle in a quartic double well. The minima 
of the well are at ±go = ±(2?7) 1 / 2 . The instantons for this problem are known, and the one 
that runs from q to —go is given by 

q(t) = —q tanh'yt, (2.7) 
p(t) = — im^qQ tanli7t sech 2 ^t. (2.8) 

Here, t is an imaginary time, and 7 = (2fc 1 /c 2 ^) 1 ^ 2 </- If we transcribe these results back into 
J x , J y , and J z , it is manifest that the spin trajectory goes continuously from one minimum 
of the classical energy to the other without any jumps. There is nothing in the physical 
character of the problem to indicate a qualitative change as H is reduced, so as long as one 
is in the field range where the gap does not oscillate (h > h*), the instantons should be 
continuous. 

In our second argument, we study the instantons for all H. Since these lie on the complex 
unit sphere, it is highly advantageous to use stereographic coordinates, 

z = tan |0e^, z = tan \9e~^ : (2.9) 

where 9 and (f) are the customary spherical polar coordinates. If a direction n has coordinates 
(9, 0), and if |n) is the spin coherent state with maximal spin projection along n, then up to 
normalization and phase, |n) is identical to \z) = e zJ ~\ J, J) where J 2 | J, J) = J (J + 1)| J, J) 
and J Z \J, J) = J\J, J). A point on the complex unit sphere is parametrized by complex 9 and 
<f), which means that z and z need not be complex conjugates; rather, they are independent 
complex variables. Either coordinate system shows that the complex unit sphere is a four- 
dimensional manifold, but the (z, z) system handles coordinate singularities better. In terms 
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of z and z, the cartesian components of J are given by 

J 



[z + z, —i(z — z), 1 — zz) 



(2.10) 



1 + zz 

The south pole is the point z = z = oo; if z — > oo with z finite, J — > J(z~ l , — iz -1 , —1); if 
z — ► oo with z finite, J — > J(z _1 , zz" 1 , —1). 

We have found the instantons for this problem in Ref. [?| among other places, but we 
restate the results along with some attendant formulas for ready reference, and to enable 
readers to verify the key results for themselves. Let us define 



E(z,z) = (z\H\z)/(z\z). 

Then the Euler-Lagrange equations of motion for the instantons are given by 

^ (l + zz) 2 dE . (1 + zzfdE 



(211) 



2J dz 



2J dz 



(2.12) 



With the Hamiltonian (|1.1|) 



E(z,z) = hJ 2 



(1 - zz) 2 - X(z - z) 2 - 2h(l - z 2 z 2 ) 



1 + zz 



5^2 



(2.13) 



so that 



k\ J 



{1 + zz 



-2z(l - zz) + X(z - z)(l + z 2 ) + 2hz{\ + zz 
2z(l - zz) + X(z - z){l + z 2 ) + 2/iS(l + zz 



(2.14) 
(2.15) 



(1 + zz) 

The minima of E are at z = z = ±z where z = [(l-/i)/(l + /i)] 1/2 . There are two solutions 
to the equations of motion that run from zq to — zq. The first is given by 



z«(f) 
z«(t) 



-Zq tanh r, 

\/Atanhr + y/l - h 2 
v / A + v 7 ! - /i 2 tanhr' 



-zq 



(2.16) 
(2.17) 



while the second is given by 



z^{t) 



\/AtanhT — \f \ — h 2 
\/A — \/l — /i 2 tanhr' 
Zotanh r. 



(2.18) 
(2.19) 
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We have defined r = ut/2, with uj = 2k l J[\(l - h 2 )} 1 / 2 . Readers can verify by direct 
substitution into the equations of motion that these are indeed solutions. Also, we have 

z V)(t) = -z^(-t), z^ 2 \t) = -z^(-t). (2.20) 

For h < h*, the two solutions can be seen to be distinct by noting that in the first z 
diverges at some t with z remaining finite, while in the other the converse happens. The 
action for these two instantons differs by a field dependent phase, which gives rise to the 
gap oscillation. 

For h > h*, on the other hand, the two solutions are physically identical. To see this, let 
us note that we can write 

z (1) = -2 tanhr, z {1) = -z tanh(r + r ), (2.21) 
z (2) = -2 tanh(r - r ), z {2) = -z tanhr, (2.22) 

where, 

i /l - h 2 

r = tanh" 1 J - . (2.23) 

If h > h*, To is real, so the two solutions are related by a simple translation in time. Such 
translations are automatically included in the sum over multiinstanton paths, so there is 
only one type of instanton. 

More importantly, the instanton is entirely continuous. We plot the components of J 
versus r in Fig. [p. As can be seen, there is no discontinuity in any component. This figure 
should be compared with Fig. 2 of Ref. |4], which pertains to the problem with the C term 
in Eq. (|1.4| ), and shows instantons with patent jumps. 

Thirdly, let us consider the problem in spherical polar coordinates, as is also done by 
CH, and in Ref. |[ CH proceed by integrating out cos#, which is feasible as the integral 
is Gaussian. The resulting action for is like that of a massive particle. If one now seeks 
instantons for this problem analogously to the particle in a quartic double well say, and 
insists on a solution 0(r) in which <p is real, one discovers that the variable r cannot be 
real, but rather must follow a Z-shaped contour in the complex plane, consisting of two 
segments parallel to the real axis, joined by a segment along the imaginary axis 0. This is 
very different, however, from a discontinuity in <fi itself. Indeed, CH seem aware of this point, 
for they say "It is clear from the shape of the potential that all three parts of the trajectory 
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FIG. 1: Components of the complex vector J(r) for the instanton solution, versus r. The x and z 
components are real, the y component is pure imaginary. 

join smoothly, because 4>(t) and 4>(t) coincide at the joints. " We agree with this. And, it is 
clear that if one finds the dominant path for cos6 l (r), i.e., the value of cos# for which the 
original Gaussian factor is maximal, one will discover that cos#(t) is also continuous. To 
eliminate any doubt, we show this by examining the explicit solution for the instanton from 
Ref. ^ In our present notation, this reads 




COS( 



cos 9 



_. (/i 2 + A-l) 1/2 tanh2r 
~*(l-fc 2 -Atanh 2 2r)V2 

h{} -h 2 - Atanh 2 2r) - (1 - /i 2 )A 1 / 2 (/i 2 + A - l) 1 / 2 sech2r 
1 - h 2 - A + h 2 \ sech 2 2r ' 



(2.24) 
(2.25) 



From these one can construct the combinations tan |0e ^ to obtain z and z. After a certain 



amount of algebra, we obtain 



z(t) 

m 



-zo 



tanh r ± a 
1 ± a tanh r ' 



where 



If we now note that 



a 



{h 2 + A - I) 1 ' 2 - A 1 / 2 

(i - h 2 yi 2 

tanh -1 a = — To/2, 



(2.26) 



(2.27) 



(2.28) 
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we see that Eq. ( |2.26| ) is the same as Eqs. ( |2.16|) and ( |2.17|) up to a shift in r. Again, we 



see that there is no discontinuity in the instanton. 

It is a separate issue as to what contour what chooses for r. If, instead of the real line, 
one lets r run on the Z-shaped contour in Ref. [], z will run from z to — z along some 
contour other than the real line, while z will run along some other contour, such that 

[This can be verified directly from Eq. ( |2.26|) , and it is a consequence of energy conservation 
along the instanton.] Since the kinetic term in the action can be written as a contour integral 
over z, 



S = J 



1 



dz 

z{z) ~ z Tz 



, dz (2.30) 

1 + zz 

it follows that any deformation of the z-contour that does not cross any singularities of z(z) 
is allowed. This is indeed the case for us, as for h > h*, the pole at z = + h) lies 

outside the segment [—z , zo] on the real line. 

We conclude by clarifying the misconception in Ref. [| The discussion above shows that 
there is only one instanton when h > h*, whereas, in the earlier work it is implicit that there 
are two, with equal, real Euclidean actions 0. However, the central thesis of that work, 
and the calculation of the action itself, are completely unaffected. 

One further point is that at exactly at h = h*, standard instanton calculations of the 
Hamiltonian (|1 . 1| ) cease to be valid, because two distinct instantons merge into one. This 
is rather like the phenomenon of merging critical points in the steepest descent approxima- 
tion for ordinary integrals |[L0|| . We expect existing formulas for the tunnel splitting to be 
inaccurate in a narrow window of fields near h = h*. The gap can undoubtedly be bridged 
via a uniform asymptotic approximation, which could be developed using the discrete phase 
integral method as in Ref. |TT[ In the methodology of that paper, what happens as we cross 
the point h = h* from below is that the width of the unusual forbidden region, in which the 
wavefunctions are exponentially damped with oscillations, shrinks to zero. In this case, one 
should join the wavefunctions on the two sides of this region using one quadratic turning 
point connection formula as opposed to two linear turning point formulas. We leave this as 
an exercise for the reader. 

In summary, we have shown that that the high-field instantons in the Feg problem are 
continuous, and thus qualititatively distinct from the discontinuous instantons induced by 
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the fourth order anisotropy in Eq. (|1.4|) . 
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